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ABSTRACT 

We investigate self-consistent particle acceleration near a pulsar polar cap (PC) 
by the electrostatic field due to the effect of inertial frame dragging. Test particles 
gain energy from the electric field parallel to the open magnetic field lines and lose 
energy by both curvature radiation (CR) and resonant and non-resonant inverse 
Compton scattering (ICS) with soft thermal X-rays from the neutron star (NS) surface. 
Gamma-rays radiated by electrons accelerated from the stellar surface produce pairs 
in the strong magnetic field, which screen the electric field beyond a pair formation 
front (PFF). Some of the created positrons can be accelerated back toward the surface 
and produce 7-rays and pairs that create another PFF above the surface. We find that 
ICS photons control PFF formation near the surface, but due to the different angles 
at which the electron and positron scatter the soft photons, positron initiated cascades 
develop above the surface and screen the accelerating electric field. Stable acceleration 
from the NS surface is therefore not possible in the presence of dominant ICS energy 
losses. However, we find that stable acceleration zones may occur at some distance 
above the surface, where CR dominates the electron and positron energy losses, and 
there is up-down symmetry between the electron and positron PFFs. We examine the 
dependence of CR-controlled acceleration zone voltage, width and height above the 
surface on parameters of the pulsar and its soft X-ray emission. For most pulsars, we 
find that acceleration will start at a height of 0.5 - 1 stellar radii above the NS surface. 
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1. INTRODUCTION 

The theory of particle acceleration in pulsar magnetospheres has been under development for 
almost three decades. Although it was well known that rotating magnetic dipoles would induce 
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electric fields in vacuum (Deutsch 1955), it took several years after the discovery of pulsars to 
realize that the huge vacuum fields could not in practice be available for particle acceleration. The 
electric field parallel to the magnetic field is at least partly screened by particles supplied from 
the stellar surface (Goldreich h Julian 1969) or by electron-positron avalanches (Sturrock 1971). 
The true accelerating voltage of a pulsar must be determined by departures from the corotation, 
or Goldreich-Julian charge density, that could completely screen the parallel electric field. Several 
types of models have studied pulsar acceleration due to charge deficits at different locations in 
the magnetosphere. Polar cap (PC) models consider the formation of a parallel electric field in 
the open field region near the magnetic poles, while outer gap models consider accceleration in 
the outer magnetosphere, near the null charge surface (see Mestel 1998 for the most recent and 
comprehensive review of pulsar electrodynamics). Ruderman & Sutherland (1975; hereafter RS75) 
introduced a PC model invoking a vacuum gap due to the trapping of ions in the neutron star 
(NS) crust. The calculations by Jones (1985, 1986), and Neuhauser et al. (1986, 1987) seem to 
favor a low value for the work function (at least a factor of 10 less than it was thought earlier) 
in the NS surface with a strong magnetic field. The important implication of this study is that 
the possibility of free ejection of charges (actually of both signs) from the NS surface can be 
now, at least theoretically, justified. In this paper, we concentrate on a space-charge limited flow 
model (implying very low work function in the NS surface) based originally on the work of Arons 
& Scharlemann (1979; hereafter AS79), who determined the electric field produced by the small 
deparature from the Goldreich-Julian charge that grows above the surface due to the geometry of 
the open dipole field. The electric field accelerating electrons in this model developed along only 
field lines that curved toward the rotation axis ("favorably" curved field lines), so that acceleration 
occurred over half of the PC. The parallel field is shorted-out at a height above the surface 
where the 7-rays from accelerated particles produce sufficient electron-positron pairs in the strong 
magnetic field. The accelerating potential is thus limited by such a pair formation front (PFF). 

These initial calculations of electron-positron PFFs assumed that the primary electrons began 
accelerating at the NS surface and that curvature radiation (CR) was the only mechanism for 
providing pair-producing photons (Arons 1983; hereafter A83). In recent years, it has become clear 
that inverse-Compton scattering (ICS) of soft thermal X-ray photons from the hot NS surface by 
the primary electrons is also an important mechanism above the PC. As well as being a significant 
energy loss (Kardashev et al. 1984, Xia et al. 1985, Daugherty & Harding 1989, Sturner 1995; 
hereafter S95) and radiation (Sturner & Dermer 1994) mechanism, ICS can also provide photons 
capable of producing pairs. Pulsed X-rays have been detected from a number of pulsars which are 
consistent with blackbody spectra at temperatures around 10 5 - 10 6 K (Ogelman 1995). Zhang & 
Qiao (1996) and Zhang et al. (1997; hereafter ZQLH97) first explored the effect of the pairs from 
inverse-Compton photons on the acceleration in a Ruderman-Sutherland type model. They found 
that ICS photons may produce a PFF sooner (at a lower altitude) than the CR photons would 
from the same accelerating electrons. In fact in this case, the electrons will stop accelerating 
before they can emit significant CR. The standard models of PC acceleration thus need substantial 
revision. 
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Another effect which has never been included in PC acceleration models is the formation 
of a lower PFF due to the positrons that are turned around and accelerated downward from 
the electron PFF. Although the number of positrons which are accelerated downward is small 
compared to the number of primary electrons and even to the charge deficit near the upper PFF, 
the multiplicity of the downward cascades is quite large (as we will discuss in Section 3.1). Thus, 
the amount of charge produced by only a small number of downward moving positrons may be 
sufficient to establish a second PFF. Although downward going cascades have been discussed 
in previous papers (see e.g. AS79), their effect on the acceleration of primaries has not been 
investigated. Daugherty & Harding (1996, hereafter DH96) qualitatively discussed the effects of 
pair cascades by returning positrons, their creation of pairs within the acceleration zone and the 
need for a self-consistent model of PC acceleration. 

In this paper, we present a detailed study of the acceleration of primary electrons and 
secondary (downward-moving) positrons above a pulsar PC, assuming space-charge limited flow 
(free emission) of particles from the NS surface (see Harding & Muslimov 1998, for a review). 
We include the general relativistic effect of inertial frame-dragging, which induces a much larger 
electric field than that expected in the flat spacetime and is not limited to favorably-curved field 
lines (Section |2.l| ). This is important because, as has been concluded earlier in papers by Fawley, 
Arons & Scharlemann (1977) and A83, the potential drops (derived for flat space-time) are not 
sufficient to account for oberved pulsar 7-rays. Both electrons and positrons suffer energy loss 
and emit photons from CR and ICS. The treatment of ICS of both upward and downward moving 
particles requires revisions from previous studies of only upward moving particles (Dermer 1990; 
hereafter D90; S95), which are presented in Section 2.2.2| . We then compute the location of both 



electron and positron PFFs due to one-photon pair production as a function of magnetic colatitude 
and height of the lower PFF. We also discuss in Section the fraction of positrons that are 
turned around at the upper PFF. When the electrons are assumed to accelerate from the NS 
surface, we find that ICS photons produce the PFFs (Section 3.1), in agreement with the results 
of ZQLH97. However, we also find that there is substantial difference between the scattering of 
upward-going electrons and downward- going positrons by the same thermal X-ray photons: the 
electrons scatter these photons at angles less than tt/2, while the positrons scatter the photons 
head-on. Electrons produce pairs through resonant scattering, while positrons produce pairs by 
scattering above the cyclotron resonance. The photons scattered by positrons are therefore more 
energetic and produce pairs in a shorter distance. These pairs may screen the accelerating field 
up to some distance above the surface. We find that stable, double PFFs can form only when 
CR photons produce them, i.e. at a height where CR losses overtake ICS losses. We compute 
the height of these stable PFFs as a function of pulsar period and surface field strength (Section 
3.1). One interesting result is that the acceleration voltage limited by CR-controlled PFFs is 
only a function of magnetic colatitude (i.e. geometry of the open field lines), ranging between 
~ 10 7 and 3 x 10 7 mc 2 , and is insensitive to pulsar parameters such as period and surface value 
of the magnetic field strength and even to the height of the acceleration. The stable location of 
the lower PFF depends primarily on surface magnetic field, temperature and size of the hot polar 
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cap, ranging between 0.5 and 1.0 stellar radius, but is insensitive to period. Implications of these 
results for high-energy pulsar emission are discussed in Section 0. 



2. CALCULATION OF PAIR FORMATION FRONTS 

We consider a test particle approach to the determination of PFF locations. The test particles 
in this case are electrons or positrons that gain energy through electrostatic acceleration and lose 
energy through radiative losses. Photons are created by CR and ICS and destroyed by magnetic 
pair creation. We assume that the electric field is completely screened at the point where the first 
pair is produced. This is a good assumption for several reasons. First, as discussed in more detail 
below, the electric field arises due to a small imbalance between the actual charge density and 
the local, rotation-induced, Goldreich-Julian charge density. It therefore does not require much 
additional charge to short-out this field. Second, the onset of pair cascading occurs very quickly 
(DH96), so that the number of pairs produced per primary particle increases rapidly over small 
distances. Thus, as found in A83, the width of the PFF (the screening distance of the electric 
field) is very small compared to other dimensions of the problem. Figure 1 illustrates the geometry 
of the calculation. Suppose that the PFF results from pairs produced by 7-rays of energy e m ; n 
(in units of mc 2 ) radiated by particles of energy 7 m i n - The distance of the PFF from the starting 
point of the particle acceleration at ho is then: 

So 

— min[»S f a (^min) ~t~ 5p(£min)] (1) 

where S' a (7 m i n ) is the distance required to accelerate the particle until it can radiate a photon of 
energy e m i n , and S p (e m m) is the pair attenuation length of the photon. The acceleration distance, 
S a {lmin), is determined by first integrating the equation of motion of the particle to determine its 
energy as a function of its pathlength s: 

d'y e 

where En is the electric field induced parallel to the magnetic field, j IC and 7 CR are the loss rates 



for ICS and CR. Discussion of these processes will be given in Sections |2.l| and |2,2| . The pair 
production attenuation length of photons radiated by the particle through either ICS or CR is 
then determined. This attenuation length S p (e), defined to be the path length over which the 
optical depth is unity, is given by 

r(e) = f SA£) T pp (d kB ,e)ds = l , (3) 
Jo 

where ds is the pathlength differential along the photon momentum vector k , T pp is the 
attenuation coefficient for one-photon pair production and #kB is the angle between k and the 



local magnetic field direction. Computation of S p (e) will be discussed in Section 2.2. As we will 



discuss in Section 3A, when ICS causes the PFF + and PFF~, S£ is smaller than S~ , due to 
different modes of ICS. The PFF + will then not coincide with the start of the electron acceleration 
at ho. 
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2.1. Polar Cap Electrodynamics 

We use the electric field due to inertial frame dragging above the NS surface, first calculated by 
Muslimov & Tsygan (1990, 1992; hereafter MT90, 92). The regime under which the generation of 
this field occurs is actually the same as implied in the electric field computations in flat spacetime 
by AS79, that of space-charge limited flow. An electric field must be present above the NS surface 
because as charges flow along open magnetic field lines, the corotation, or Goldreich- Julian charge 
density pqj, cannot be maintained. Even though p = pqj and therefore En = at the surface, the 
curvature of the field lines causes the area of the open field region, through which the particles 
flow, to increase faster than pqj, and a charge deficit grows with distance. Thus, the En grows 
with height up to about one stellar radius above the surface, and then drops off. It is important 
that for a nearly aligned NS in flat spacetime, the space-charge density of the outflowing particles 
(electrons) proved to be almost exactly compensated by the effective Goldreich-Julian charge 
density, thus resulting in the substantial suppression of the electric field in the region of open 
field lines. General relativity causes a significant modification of the mechanism of the electric 
field induction, through the effect of dragging of inertial frames, a consequence of the distortion 
of spacetime by a rotating gravitating body. An observer near a rotating mass experiences a 
force and must corotate to remain in an inertial frame, with an angular velocity that decreases 
with distance from the rotating mass. The charge density above a NS surface must be computed 
in a local inertial frame that is rotating with respect both to the NS and to an observer at 
infinity. Thus the Goldreich-Julian charge density, which is the charge density required to make 
magnetospheric particles drift in corotation with the star, will differ from that in flat spacetime. 
This charge difference enhances En over what it would be in flat spacetime, by a factor of 50 - 100 
for a typical 1 s pulsar. The frame-dragging contribution to En that signifies the striking difference 
between the general relativistic and classical treatments of NS electrodynamics is proportional 
to cos x (see equation (T^j below) , where x is the angle between the magnetic and spin axes of 
the pulsar. [Obviously, there is also a frame-dragging contribution to En which is proportional to 
sinx, but it is comparable to that produced in a flat spacetime limit.] Particle acceleration may 
therefore occur throughout the entire open field line region, with the relative contribution of the 
frame dragging component to En being strongest for pulsars with small obliquities. 

Muslimov & Harding (1997; hereafter MH97) derived expressions for the En due to frame 
dragging in two limits: close enough to the NS surface such that z = s/R -C 9q, where 9q is the 
PC half-angle, and far from the surface such that z 9q. These expressions were derived from 
solutions to Poisson's equations assuming boundary conditions En = 0, as well as the potential 
$ = at the stellar surface and along the last open field line. We have adapted these solutions for 
use in this paper by incorporating the screening effect of an upper boundary at z c = S c /R, i.e. of 
the pair formation front, where En{z = z c ) = 0. Although this approach is not fully self-consistent, 
in Section 2.1.1 we discuss how the screening of the electric field at the upper boundary could be 
included in our calculations in a self-consistent way. In this paper we also explore the situation 
where the positrons flowing back to the PC surface initiate electron-positron cascades. Although 
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this possibility has been discussed since the very first papers on pulsar electrodynamics (see 
e.g. AS79), it has never been addressed at an appropriate quantitative level. Here we attempt 
to approach this and some other related problems from the quantitative point of view which, 
we believe, can advance our understanding of the PC electrodynamics and its relevance to the 
modeling of 7-ray emission from pulsars. 

We note that the derivation of formulae for the electric field (potential) we exploit in the 
present analysis implies that the angle #0 (magnetic colatitude of the polar field line) is small, 
which is a very good approximation for the region of open magnetic field lines in most pulsars. 
However, it is instructive to look at the general-relativistic expression for the Goldreich- Julian 
charge density that contains a contribution of order of 6 2 to the main term. This expression reads 

cosx + — -&(?/)# sin x cos (/>j- > (4) 

where f2 is the NS rotation frequency, Bq is the surface value of the magnetic field strength at the 
magnetic pole, a = (1 — e/rf) 1 ^ 2 is the red-shift function, e = r g /R, r g is the gravitational radius 
of the NS, R is the stellar radius, n = r/R is the dimensionless radial coordinate, k = el /MR 2 , 
I and M are the moment of inertia and mass of the NS, respectively, and the functions f(rj) and 
H(rj) are defined below (see equations [8], [16]). In formula (4) the functions / and H should be 
evaluated at rj = r/R and e = r g /R. In our derivation of the general-relativistic electrodynamic 
equations we neglect the contribution to the electric field produced by the flaring of the magnetic 
field lines, which is a factor of ~ 9~ 2 ~ 500 (P/0.1 s) (where P is the pulsar spin period) smaller 
than that produced by the frame dragging. This contribution results from the second-order terms 
like that in the square bracket of equation [4]. [In Section 2.4 we will refer to this equation to 
clarify yet another important issue.] The small-angle approximation may not be accurate enough 
for the millisecond pulsars, for which 9 ~ 0.3 (r / R) 1 / 2 (P/2 ms) -1 / 2 . We also note that the effect of 
the electric field generation by the frame dragging does not depend on the particular configuration 
of the stellar magnetic field, simply because while the distribution of the real space charge in the 
acceleration region is mainly determined by the flaring of the magnetic field lines, the effect of 
frame dragging on the Goldreich-Julian charge density is independent of the geometry of the polar 
field lines. This means that one cannot mimic the frame dragging effect simply by distorting the 
polar magnetic flux tube. 



2-Kcari 3 f(l) 



-H(t])6 2 



2.1.1. Reseating of Ej| with no upper PFF 

The possibility of electron-positron cascades near the stellar surface initiated by the 
backflowing positrons unavoidably implies that the bottom of the polar magnetic flux tube 
should be treated as a highly-conducting layer of an electron-positron plasma on the top of the 
PC surface, rather than as a regular NS surface made of iron or hydrogen atoms in the strong 
magnetic field. In a more extreme case (see e.g. Wang et al. 1998), which seems to be less justified 
at the moment, the whole NS surface may be covered by a thin electron-positron plasma layer. 
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Note that the presence of such a layer even on the top of the PC should substantially affect the 
electrodynamics and therefore the acceleration of charged particles in the region of open magnetic 
field lines. It will allow the free supply of charges into the acceleration region, a requirement for 
space-charge limited flow. Also, we shall discuss in Section 3, a self-consistent regime with double 
PFFs that may significantly differ from the standard one implying a single PFF above the stellar 
surface. 

In this Section and the next, we discuss the principal modifications one must introduce into 
the treatment of the acceleration region to incorporate the effects resulting from the occurrence 
of double PFFs. First of all, the appropriate distance variables in En should be rescaled by an 
"effective" NS radius -Re = R + ho (cf. Fig 1), which is now assumed to be the lower boundary of 
the acceleration, such that E»iz = 0) = 0, where z = s/Re is the dimensionless altitude above the 
effective surface. The expressions for the electric field are thus symmetric between the lower and 
upper boundaries, in the sense that they allow the proper treatment of e.g. both primary electron 
and secondary positron acceleration, which is an essential requirement for the theory. 

To clarify our analytic solutions (see Section 2.1.2) for the rescaled electric potential 
(accelerating electric field) produced between the lower and upper pair fronts, let us consider 
the expression for the electric potential in terms of the difference between the actual and 
Goldreich- Julian charge densities p — pcj. Because p and pqj are just linear combinations of sinx 
and cos x, we can write the potential $ in the simple form (at distances from the stellar surface 
greater than the PC size) as 



$ = aS(r) 



(p ~ Pgj)| x=0 cosx + gG - PGj)| x= ^ sinx 



(i-n, (5) 



where S(r) = ir[r8(r])] 2 is the cross-sectional area of the polar magnetic flux tube at the radial 
distance r and (p — pgj)\ x = and (p — pgj)\ _ 1L are the coefficients of cosx and sinx in the 



expression for (p — pgj) (see equations [ 10 1 and |TTJ). It also says that the local value of the electric 
potential is of order of the local value of the effective space charge (equal to the difference between 
the real space charge and induced Goldreich-Julian space charge) divided by the characteristic 
longitudinal length scale. The expression (]5|) does not incorporate the effect of the screening 
of the electric field at the upper boundary (PFF). This effect is included in the corresponding 



expressions we discuss later in Section 2.1.2j . The relation (|5|) gives the correct solution for the 



electric potential and longitudinal component of the electric field for altitudes greater than the 
PC size and satisfying the $ = and En = boundary conditions at the (effective) stellar surface, 
where p = pgj- This expression is very useful for illustrating the modification of the potential and 
longitudinal component of the electric field resulted from the redefining of the position of the lower 
zero-electric field boundary (set e.g. by the lower PFF produced by the backfl owing positrons). 
It explains the physical meaning of such a modification, which is simply a readjustment of the 
true charge outflow (in the regime of self-limitation) to the local value of the Goldreich-Julian 
charge density at the effective surface, and allows us to better understand the corresponding exact 
analytic solution. 
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We introduce the radial coordinate scaled by effective radius -Re, rj = 1 + z = (Re + s)/Re, 
and transverse coordinate £ = 0/6(rj), the magnetic colatitude scaled by the half-opening angle of 
the polar magnetic flux tube, 

1/2 



0(7/) = O 



at radius 77, where 



?0 



is the PC half-angle at the effective surface Re, and 

3 r 



1/2 



/(,)-. (2)' 



In 1 



»7 



+ 



77 



1 + 



277 



(6) 
(7) 

(8) 



is the correction factor for the dipole component of the magnetic flux through the magnetic 
hemisphere of radius r in a Schwarzchild metric (MH97). 
Note that in equation (5) 

Q r 3 

The expressions for p and pc,h valid for any radius, can be written as (see e.g. MH97) 



p = -o-(r) 
-a(r) 



(1 — rjln) cos x + ~ &oH(l)^ sin X cos < 



(1 - — )cosx + - 6 (r)) H(rj)£ sin x cos i 

7? Z 



In these expressions 



i? 3 



(10) 
(11) 

(12) 



/(?7*) aS{r) 

where 77* = R/Re (0.5 < 77* < 1, as it will be shown below). 

Thus, after substituting the above expressions for p and /Oqj into equation (5) we get the 
following formulae for the rescaled electric potential 



$ = 



^ o 0qV* \^l K (l - ^3) cosx+ 

\ [0(r))H(n) - 6 H(1)] e sin x cos </>} (1 - £ 2 ), 



(13) 



and the accelerating component of the electric field 



2^,2 



->ri* cos x + ^-9(ri)H(r))5(r))t sin x cos 



(i-a, 



3 9 k 

.2 77 

where $0 = Bo(£lR/c)Rf(l)/ /{r}*), and i?o = $o/R- The quantities e and k in all expressions 
derived in this Section for the rescaled radial coordinate now read 



(14) 



e = 



2GM 



Re c 



2 • 



k = 



el 
MR 2 ' 



(15) 
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Finally, the functions 



TT . e k ( 3e k 



d_ 

drj 



\n[H(r,)e( V )] 



(16) 
(17) 



are both evaluated at 77 = r / R-&. 



2.1.2. Screened E11 with upper PFF 



We now include the effect of the electric field screening at the upper PFF, at 
77 = f] c = (R + H c )/Re, by obtaining the solution to Poisson's equation with the upper 
boundary condition, #11(77 = r/ c ) = 0. In the limit where z c <C 1 and z < z c the accelerating 
component of the electric field reads 



#11 ~ -#o^o(l -t) 1/2 v. 
" 00 



2 1 2 



Li=l 



cosx+ 



.1=1 



sin x cos < 



(18) 



where, 



Ai- -K 



kfJi(h) 



^i(z,7i), B i = -e H (1)5(1) 



16 



KMi), (19) 



T(z,i) = - [01(777- l)e lz + 02(777 + l)e 72 +ai(l -7) -a 2 (l +7)]/(«i + a 2 ) 
ai = (7% + l)e" 7Zc - 7-1, a 2 = 7 - 1 - ( 7 t? c - l)e^ c , 



and 



ki 



li 



G (1 -e)V2' 



and 7i 



(20) 
(21) 

(22) 



where fcj and ki are the positive roots of the Bessel functions Jo and J±, respectively. In expressions 
(20), (21) 7 = 7j or ji should be used. In the limit, where z = 77 — 1 » 9q an d z c = 7/ c — 1 » 6>o 
(cf. Appendix A, equation [A4]), 



^11 - - E M^ri; 
+l6( V )H{f,)9{r,) 



V_\ 3 >^ 8J (kiQ {Vc)(Vc _ v) 



y, 16Jl(fcjg) e _^. (r)c)fe _ ??) 



VcJ ~i kfJ^ki) 
_ Vc O(Vc)S(Vc)H(r] c 
' 77 ^(77)0(77)^(77) 

sin x cos > 



cosx 



(23) 
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where now, 

li= 0(7,^(1-6/170)1/2' and ^=0(^(1-^)1/2 ( 24 ) 

Simple analytic expressions can be derived in several limiting cases some of which are summarized 
in Appendix A (for the radial coordinate scaled by the true stellar radius). The solutions ( |i"8|) and 
(|23| ) thus incorporate both the effect of the rescaled lower boundary and screening of the electric 
field at the upper boundary, and we will use them in our numerical calculations described later on 
in Section (3.1). 



2.2. Radiation Production and Losses 



In this section, we describe our treatment of the radiation processes that affect the polar 
cap particle acceleration. This includes energy losses due to CR and ICS, and pair production 
by the photons from these processes. Although the particles are radiating a full spectrum of 
photons, we are concerned here only with the pair producing photons. In the interest of making 
our numerical code as efficient as possible (but with some loss of accuracy), we do not model the 
entire radiation spectrum of each process, but compute pair production attenuation lengths for a 
single representative photon energy at each step along the particle path. 



2.2.1. Curvature radiation 

The CR loss rate for particles of charge e moving along a magnetic field with radius of 
curvature p c is 

3 e 2 c 

-^ = W^ 4 - (25) 
The radius of curvature of the magnetic field in a Schwarzschild metric is 

™ d GW = tviZ - £' < 26 > 

in the small angle limit, where x = r g /r, r g = 2GM/c 2 is the gravitational radius, and the function 
f(x) has been defined in equation (||). Formula ( p6| ) gives a slightly larger radius of curvature than 
the flat-space formula does (by a factor of G, about 25 — 30%, for the NS parameters we assume in 
this paper). This can be easily understood, since the strong gravity of the NS tends to increase the 
flaring of the polar field lines, and therefore the last open field lines (i.e. those reaching the light 
cylinder at their "turning-points") should emanate from slightly smaller magnetic colatitudes. In 
other words, the polar field lines get slightly straighten out, and the effective PC radius slightly 
decreases (by the same ~ 25 — 30%, see also Gonthier &: Harding 1994). 
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The spectrum of CR photons from a particle with energy 7 is a power law at low frequencies, 
with an exponential decline at high frequencies: 




™ C r0o<* ; (27) 



where e CI = (3/2)(7j/mc)7 3 /p c is the critical frequency. Since the one-photon pair attenuation 
coefficient (equation |48| ) increases sharply with energy as exp (— 8/3eB' sin#kB), the photons 
from the exponential tail of the curvature spectrum will produce the bulk of the pairs. Here, and 
in subsequent Sections, the local dipole magnetic field strength is 



B' 







) 3 


\ 7(1)1 


CO 


{re, 




lf(v*)\ 



(28) 



where B CT = 4.413 x 10 13 G is the critical field strength. The overlap of the curvature spectrum 
and the pair attenuation coefficient will fall in a narrow frequency band, which we can approximate 
in a steepest descents analysis as a Gaussian of width, Ae p = (3B'e p sin^kB) 1 ^ 2 , with a mean 
energy e p = (8£ cr /3B' sin O^b) 1 ^ 2 . Photons will pair produce roughly when e p B' smO^B ~ 0.2 for 
B' < 0.1 and as soon as threshold £ p sin#kB = 2, is reached, for B' > 0.1 (Daugherty & Harding 
1983; hereafter DH83). We therefore have 

e p = 13.3£ cr , Ae p = 0.273e p , B' < 0.1 

£p = 4e cv /3B', Ae p = 0.87B' 1 / 2 e p , B' > 0.1 1 ' 

At each step along the particle acceleration path, pair attenuation lengths are computed for up to 
five CR photons with energies evenly distributed between e p and e p + 4Ae p . The first finite value 
(if any) of the pair attenuation length is taken as S p for that step. 



2.2.2. Inverse Compton scattering 

As the particles are accelerated, they may scatter soft photons from the hot NS surface. 
Due to the strong magnetic field, the ICS cross section is resonant at frequencies where electrons 
may be excited to higher Landau states and is strongly suppressed (for extraordinary mode 
photons) below the fundamental. Although the full QED scattering cross section for this process 
has been calculated (e.g. Daugherty & Harding 1986, Bussard et al. 1986), it is quite unwieldy 
for numerical calculations. However, the magnetic scattering cross section in the Thomson limit 
(Canuto et al. 1971) is much simpler and has been used in place of the full QED cros section in 
almost all astrophysical applications. It includes only the resonance at the cyclotron fundamental 
(i.e. scattering in the ground state with no excitation) and approaches the Thomson cross section 
at frequencies above the cyclotron energy. 

In this paper, we follow the treatment of S95, who derived simple expressions for electron 
energy loss rates for magnetic Compton scattering in the Thomson limit, expanding the treatment 
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of D90 to include the Klein-Nishina suppression above the resonance. Dermer's approach was to 
divide the Compton scattering loss rates into three components, each due to the different behavior 
of the cross section in the regions below the resonance (which he called "angular" scattering 
due to the strong angle dependence), in the resonance ("resonant" scattering) and above the 
resonance ( "non-resonant" scattering) . He assumed that the resonant part of the cross section 
was a (5-function at the cyclotron energy, and that the soft photon source was a blackbody of 
temperature Tb radiated by a hot PC of radius Rt- He derived the scattering loss rates for 
electrons moving away from the NS surface, so that the soft photon density decreases with height, 
due to a dilution factor. We have expanded this treatment by computing the associated loss rates 
for particles (in our case, positrons, although the sign of charge makes no difference for Compton 
scattering) moving toward the NS surface. In this case, the particles approach the soft photons 
"head-on", scattering them to higher energies than the downward moving particles. 

The Compton scattering energy loss rate for a particle of energy 7 moving through a photon 
field, n p h(e), that is uniform between angles /i_ < fj, < /x+, from equation (12) of D90, is 

C [["/(l-Pl-i-)}' 1 ft*+ fl f£'s,max ( da' 

~^ = ]^zJo den ph (e)J^ d^l-^J ^J^ 

(30) 

where e' s max = 75(1 — (3/j,) is the maximum scattered energy in the particle rest frame. Here, 
primes are used to denote quantities in the particle rest frame and the subscript s denotes scattered 
quantities. Thus, a' is the rest frame scattering cross section. Based on these expressions, the loss 
rates for the three parts of the cross section, 7 IC = 7 ang + 7 res + 7 KN , used in equation ( |30|) can be 
written 

- 7ang = 46.1 T 6 4 (1 - fie) / ang s -1 , (31) 



7res = 4.9 x 10 11 



T e Bf 2 



/res S , (32) 

PI 

- 7kn = 3.7 x 10 11 (Jj^j £ + de' [s + e /e' - 7] /kn s" 1 , (33) 

where Tq = T^/IO 6 K is the NS surface temperature, B12 is the local dipole magnetic field strength 
in units of 10 12 G. equation (|33|) assumes a J-function at energy eq = 2.70T, where = kT/mc 2 , 
for the distribution of soft thermal photons, and 

, _ 7^2 _ _ (//+ + //_) In [(1 -/?//_)/(! -/?/,+)] 

/ans 7 2 "1 30i+ 2/3 7 2 + /? 3 7 4 (/x + - M-) ' 1 j 



/r, 



^ ) In 



1 \ 1 



+ r 



{[w + ln(l — e 



w_ ln(l — e 



Li 2 (e- W -)}} 



(35) 
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/KN 



2(e' 2 - 4e' - 3) 2^ 



2(1 + e')(e' -2e'-l) 



4e' 2 + 6e' + 3 
(l + 2e') 3 
V 2 - 2e' - 



e'(l + 2e') 



A2 



ln(l + 2e')--. 



(36) 



In the above expressions, w± = e B /[#x7(l — /3/"±)]> £± = £q7(1 — Ph±) and £ B = B' is the local 
cyclotron energy in units of mc 2 . Equations (31), ( |32| ) and ( |33"| ) are identical to those given by 
S95, but the expressions for f an g, / res and /kn have been generalized to allow for scattering 
by downward moving particles. Specifically, the corresponding equations given by S95 have 
specialized to the case of upward moving particles (//_ = /i C) /^+ = 1) see equation ([}?]) below). 



equation (|34|) is the same as D90's equation (16), which is general enough to compute the loss 



rate of either upward or downward moving particles. Equation (35) includes, to first order, the 



rest-frame particle recoil e' s ~ e'[l — £'(n' s — //) /2] and is a major improvement over previous 
expressions for resonant Compton scattering that assume no recoil (e' s = e'). The first term in 
equation ( |35| ) gives D90's equation (54), and the remaining terms result from recoil, where Lii 
is the Dilogarithm function. The treatment of recoil in resonant scattering is necessary in this 
calculation to compute accurate energies for photons scattered by the electrons and positrons. 

Given the generality of the above expressions, the only difference in computing scattering loss 
rates for upward-moving electrons and downward-moving positrons lies in the integration limits 
over polar angle \x in the laboratory frame. For the case of semi-isotropic blackbody radiation at 
height h = s + ho above a hot PC of radius R^, this translates into the following values for 



where 



electrons: 



positrons: 



fj, c = cos 6 C 



At- = He 
M+ = 1 
M- = -l 

= -He 

h 



h 2 + Rl 



(when i?T "C R) gives the dilution factor that is present in the soft photon density, 

1 - He 



where 



^ph(e) = n hh (e) 
«bb(e) 



8tt 



\ 3 C [exp (-e/9 T ) - 1] 



(37) 
(38) 

(39) 

(40) 
(41) 



is a blackbody spectrum of temperature #t and A is the electron Compton wavelength. 

To determine the energy of the scattered photons that may produce pairs, the scattering 
kinematics of the electron and positron must be treated accurately. For non-resonant scattering in 
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the Thomson limit, i.e. when Ajeq < 1, we can assume there is no recoil in the particle rest frame, 
and the maximum scattered energies in the laboratory frame are 

e s = 27 2 (1 - fi c )e electrons, 
e s = 47 2 £o positrons. 

Since the electrons and the soft photons are both moving away from the NS surface, e s depends on 
fj, c , which decreases with height. However, the positrons scatter the soft photons "head-on" with a 
maximum angle of 7r, and e s is independent of height. When 4769 > 1, or for resonant scattering, 
we must include recoil in the rest frame to compute the scattered energies. The average scattered 
energy in the laboratory frame is computed from the scattering rates, which for 7-es 

(equation 

p2|) and 7 KN (equation [ 33 1 ) explicitly include recoil, and the loss rates, so that 



(e.) « JT- (43) 

for each component, where 

Ric = / den ph (e) / d» / d// s / de' s — — . (44) 

— yU_ Jo J-l Jo \dlJL' s d£' s J 

(e s ) is thus the weighted average of (e s — e), where (e s 3> e) for scattering by relativistic particles. 
Since photons scattered by the "angular" part of the cross section below the resonance do not 
have energies above pair threshold, for our calculation we need only R ms and -Rkn- Using equation 



(44) and following Dermer's (D90) method of computing 7 res , we have 



^ = 2.18xlO-(^j ln (^_Ij s-K (45) 

For 7 > 1, the resonant scattered photon energy, e v s es ~ jB' is independent of incident angle, and 
is thus the same for upward and downward moving particles. However, the energy of the incident 
photon in the particle rest frame, = £07(1 — /3/x±) will not be the same, so the upward and 
downward moving particles will not necessarily both undergo resonant scattering. Using equation 



(44) and following Stumer's (S95) method of computing 7kn> we have 



R KN = 3.7 x 10 11 ) / de'r KN s" 1 , (46) 



£' 



- + n \ o.n 2 + " 7T2 " ln (1 + 2e )• ( 47 ) 



where 

4 , 2e'{l +e') (e' 2 -2 e' -2 
' (l + 2e') 2 + e 

Our numerical calculations show that (e^ ) approaches 7 for 4760 3> 1, as expected. Thus, 
(e^ N ) > (ef 8 ) for B' < 1. 
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2.3. Pair Production 

The method we use to compute the electron-positron pair attenuation length of the photons, 
S p (e), has been described in detail in Harding et al. (1997). Using equation @, S p (e) is computed 
by integrating the pair production attenuation coefficient of the photon along its path through 
the dipole field. The photon is assumed to pair produce at the point where t(9,e) = 1, and 
S p (e) is then set to that path length. The two main inputs needed are the energy of the photon 
and its angle to the magnetic field, #kB- The energies of the radiated photons were given in 



Sections 2.2.1 and 2.2.2. For the angle of the radiated photons at the emission point, we assume 
that 6*kB = 2/7 for Compton scattered photons (D90) and 6*kB = for CR photons. Due to the 
curvature of the field lines, #kB will grow as the photon propagates, roughly as sin6>kB ~ s/p c . To 
evaluate the pair production attenuation coefficient at each point along the path of the photon, we 
Lorentz-transform the photon energy and local magnetic field to the frame in which the photon 
propagates perpendicular to the local field. This is the center-of momentum frame for the created 
pair, where the attenuation has its simplest form and the photon energy is e CM = esin#kB- The 
one-photon pair attenuation coefficient is considered in two regimes. For B < 0.1 B cr , photon pair 
produce far above threshold, where the asymptotic expression in the limit of large numbers of 
kinematically available pair Landau states (Tsai & Erber 1974, DH83) can be used: 

1 a 

2 A 



2fP = ~ B ' A \U (X), (48) 



(0.31,0.15) exp(-^) X «l 
A, U (\)H ' (49) 

(0.72, 0.48) x~ 1/3 X»l 

where x = / £ cm/^^' 1 a ^ s the fine-structure constant, B' is the dipole field strength at point s 
along the photon path (see equation |p8|j). When B > 0.1 -B cr , pair production will occur near 
threshold, where the above expression is not accurate. We thus include the factor / = 1 + 0.42e~^' 7 
in x, introduced by DH83, as an approximation to the near-threshold attenuation coefficient. In 
this paper we compute the attenuation length averaged over photon polarization. 



2.4. Backfiowing Positrons, Double Pair-fronts and Polar Cap Heating 

Soon after the concept of electron-positron pair production had been introduced into pulsar 
models (Sturrock 1971), it became clear that the precipitation of positrons (or, in some models, 
electrons) onto the stellar surface may have some important effects. For example, Sturrock has 
noticed that (see p. 531) "if the positrons were all returned to the surface, the resulting space 
charge would reverse the sign of the electric field at the surface, cutting off the flux of primary 
electrons". Then he added that "one possibility is that the configuration so adjusts itself that 
the primary flux is never quite cut off, so that the flow is steady, and another possibility is that 
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the flow is oscillatory" . Thus, he drew attention to the existence of an intrinsic feedback between 
the pair formation and the electric field accelerating primary particles. Later on, several authors 
(see e.g. RS75, Jones 1978, AS79, and Arons 1981) have discussed similar ideas focusing on the 
possibility of PC heating by the flow of energetic positrons (electrons). DH96 discussed the effect 
of downward moving positrons on the particle acceleration and the possibility of a lower PFF. 
The most detailed study has been performed by AS79 who calculated the fraction of returning 
positrons, structure of the PFF with positron trapping near but below the lower boundary of 
the PFF, and briefly discussed the positron-initiated cascade very near the stellar surface. They 
estimated the fraction of returning positrons as ~ 9q (<C 1, where 6o is the half-opening angle 
of the polar flux tube at the stellar surface) of the Goldreich-Julian charge density. The main 
argument underlying this estimate is that the contribution of the second term in the square 
bracket in equation (4) increases with altitude (due to the flaring of the magnetic field lines 
emanating from the PC), and that (see AS79, p. 867) "the additional negative charge density 
needed to achieve En = 0, over and above that of the electron beam, is therefore a fraction ~ 6q 
of the negative charge density already present in the electron stream". And, "as the pairs begin 
to form near but below some height (at which the electric field shuts off), the residual Ej| of 
the diode causes the secondary electron-positron plasma to polarize and shield these "external" 
charges. The formation of this polarization charge requires a dynamical response which, under 
most conditions, leads to the formation of the downward directed positron stream whose flux 
at the stellar surface is F + <C -F-". However, the main problem with this reasoning is that if 
we assume the zero-electric field boundary condition at the stellar surface, then "the additional 
negative charge density" is nothing else but the imbalance between the real and Goldreich-Julian 
charge densities that produces the electric field above the stellar surface in the first place. This 
imbalance (or charge deficit) increases with altitude and reaches some maximum value. Thus, if 
the backflowing positrons shield these "external" charges at some height, then they will also screen 
the electric field all the way down to the stellar surface, and we should unavoidably come up with 
the situation envisaged by Sturrock, where the returning positrons cut off the flux of primary 
electrons, thus resulting in the oscillatory regime. 

In this section, we explore the possibility of occurrence of a steady regime of positron backflow 
that is not capable of disrupting the current of primary electrons. We also discuss the efficiency of 
PC heating produced by the precipitating relativistic positrons and high-energy quanta. We also 
discuss the accounting of the feedback between the returning positron flux and electric potential 
distribution in the acceleration region in the self-limiting regime of particle flow. The primary 
(electron) traversing the region of the electric potential drop emits a 7— ray photon that produces 
an electron-positron pair at the altitude /i pa j r above the stellar surface. The electron-positron pairs 
keep moving forward, with some positrons of these pairs decelerating and eventually reversing 
the direction of their motion at some altitude h c . It is very likely that the same effect tends to 
occur with the positrons from both additional photons emitted by the same electron or from a 
second generation of pairs. An accurate determination of the fraction of backflowing positrons 
would require a self-consistent calculation of the screening of E» by the cascade pairs and the 
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modification of the pair energy and spatial distribution by the screened En . Short of carrying 
out such a detailed calculation here, we shall derive an upper limit on the fraction of returning 
positrons. However, since the cascades develop very fast and the number density of particles 
increases in an avalanche, the altitudes at which the positrons from the higher-order cascades get 
decelerated should be very close to but slightly above h c . The deceleration of positrons results in 
a charge separation in the initially quasi-neutral electron-positron beam, thus effectively shielding 
the electric field where the bulk of the pair front is produced. We can justifiably assume that 
En = at h » h c . The backflowing positrons slightly suppress the voltage all the way down 
to the bottom of the polar magnetic flux tube, so that the electric field En (and also potential 
<£) vanishes at a height greater than the start-off height of the primary electrons. This occurs 
because the flux of backflowing positrons is equivalent to a corresponding enhancement of the 
total electron current and therefore of the maximum number of electrons per second to be ejected 
into the acceleration region. Thus, in compliance with the zero-electric field boundary condition, 
the ejection radius should fix itself at the larger value corresponding to a higher effective "rate of 
supply of Goldreich- Julian charge" (see also expressions [4] and |5l]] proving this statement). 



Let us introduce the Goldreich-Julian current (associated with the corresponding charge 
density) or the rate of ejection of charges per element of the solid angle into the region of open 
magnetic field lines 

ac\p GJ \(r6) 2 , (50) 



where (r9) 2 = (QR/c)R 2 rj 3 {; 2 / f(rj), and dQ^ = £d£d(f) is an element of the solid angle in the PC 
region. In the rest of this Section we present our estimate of the maximum power carried by the 
backflowing positrons. Using the general-relativistic expression for p G j (equation |11|| ) we can 
write 
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1 ~ ~T I cos X + 7;H(r] E )6(r] E )^smxcos< 
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(51) 



where 1 = (nR/c) 2 B Rc/ f(l) and t/ e = R E /R = (R + h )/R = 1 + h /R > r). 



The zero-electric field boundary condition at rj = 1 requires the magnitude of the primary 
electron current flowing from the PC surface to be equal to {dX^j / dQ,^)\ . This expression 
implies (compare with expression [4]) that an enhanced electron current is ejected from the 
effective radius in the regime of self-limitation or, equivalently, allows for the additional current of 
positrons flowing from the upper PFF downward. It is important that in the steady-state regime 
the total current of actual charges (electrons and positrons) remains constant along the field lines, 
and is fixed by the Goldreich-Julian current at the zero-electric field boundary (i.e., in our case, at 
the effective radius Re)- The electric field above Re (in the regime of self-limitation) is produced 
by the imbalance between the actual current and the local value of the Goldreich-Julian current. 
At distances greater than the effective radius Re for the Goldreich-Julian current we can write 



dl GJ (r] > t/e) _ 2o 
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(52) 
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Let us consider the situation when the fraction of positrons returning from the upper PFF 
is not only sufficient to produce a lower PFF, but is also big enough to affect the electric field 
accelerating the primary electrons. For example, the backflowing positrons tend to reduce the 
accelerating electric field above the effective surface of radius Re, simply because the positron 
current is equivalent to the enhanced electron current which unavoidably requires the electrons to 
be ejected from the greater heights above the effective surface. Thus, the returning positrons may 
eventually screen the accelerating electric field and trigger an oscillating regime with alternating 
phases of developing and collapsing of double pair fronts. As has been mentioned above, the 
oscillatory regime akin to that we discuss here was anticipated in his classical paper by Sturrock 
(1971), even though he did not discuss the lower PFF. The occurrence of the oscillatory regime 
depends on the fraction of backflowing positrons, which is determined by the particle kinematics 
within the upper PFF and the penetration depth of the electrostatic field into the electron-positron 
plasma cloud. A quantitative description of this process would be possible through a detailed 
analysis of the electron-positron cascades and their feedback on the electrodynamics of the pair 
formation region. We can now estimate the maximum positron current needed to screen the 
electric field accelerating the primary electrons. For this purpose we assume that the bulk of the 
lower pair front sets up at the effective radius Re, and that the upper pair front establishes at the 
radial distance r > Re- Then, for the maximum positron current we can write, using equations 
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R/Re- The maximum value of the electric potential can be estimated as (cf. equation 
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Now we can derive the maximum total power put into the backflowing positrons (cf. MH97, 
equations [76]- [78]): 
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is the spin-down luminosity of a pulsar. A similar estimate for the maximum efficiency of 
acceleration of primary electrons yields 
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which amounts to 10 % (for a small obliquity and 77* ~ 1). 

For relatively small obliquities (or small magnetic polar angles and non-orthogonal rotator) 
and typical pulsar spin periods of 0.1 - 1 s the first term in equation ( |56| ) dominates, and we get 

{A + } max - l* 2 vl (59) 

For a 1.4 solar mass NS and for a broad range of realistic equations of state of dense matter 
(see e.g. Lorenz, Ravenhall, & Pethick 1993 and Ravenhall & Pethick 1994 for the calculations of 
the NS moment of inertia for various equations of state) I /(MR 2 ) ~ (0.2 — 0.25)(1 — r g /R)~ l . 
Thus, for the NS of 1.4 solar mass and 8-10 km radius (which is consistent with the most realistic 
stellar models) we can estimate k = (r g / R)(I /MR 2 ) w 0.15 - 0.27. Given ry* -0.5-1 (see Section 
3.1), equation (|59|) yields 

{A + } max -(3-10- 4 -2-10- 2 )(^) 2 . (60) 

This estimate combined with the recently observed X-ray luminosities (that include both pulsed 
and non-pulsed components and imply isotropic emission) of pulsars (see Becker & Trimmer 
1997) may have rather interesting implications. If the X-ray fluxes in some of these pulsars 
are dominated by the photons from the heated (e.g. by the backflowing positrons) PC, then 
luminosities higher than given by equation (|60|) would indicate that these pulsars operate in the 
oscillatory regime discussed above. According to the estimates by Becker & Triimper the pulsed 
X-ray luminosities (for the case of isotropic emission), L£, for e.g. Crab, Vela, Geminga, and PSR 
0656+14 are, respectively, 1.6 -10~ 3 , 0.7- 10~ 5 , 1.3 • 10" 4 , and 3.7 • 10~ 3 L sd (where L sd = IQ\(l\, ft 
and ft are the pulsar spin frequency and its time derivative, respectively). Although some of these 
pulsars have X-ray luminosities that lie above {-^ e +} m ax> ^ ne X-ray emission may not be isotropic 
(see Zavlin et al. 1995). If we assume that at least for Vela, Geminga, and PSR 0656+14 (see 
e.g. Harding & Muslimov 1998 for the modeling of the soft X-ray and 7— ray emission for these 
pulsars) the X-ray emission is beamed into a solid angle of ~ 1 steradian, then for these pulsars 
we can estimate that {AP} an i s = {^xlanis/^sd ~ 6 • 10~ 7 , 10~ 5 , and 3 • 10~ 4 , respectively. The 
corresponding PC temperature for these pulsars can be estimated as 

.^(o,-,.^^^.)-^) 1 ' 2 ^) 371 ^)" 374 k. m 

Thus, the estimated value of {AP} an i s ~ (3 • 10 -5 — 10~ 2 ) {A+} max may indicate that the 
backflowing positrons precipitate onto the effective area smaller than that of the standard PC 
(e.g. the returning positrons focus toward the magnetic axis as discussed below), and/or that 
the fraction of the returning positrons is well below the maximum possible one. Then the latter 
would support the quasi-steady state (MH97) rather than the oscillatory (Sturrock 1971) regime 
of pulsar operation. 

The energetics of the PC heating is mainly determined by the energetics of the backflowing 
(primary) positrons. Their energy is eventually redistributed between the high-energy photons 
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they emit and electron-positron pairs created by photons in the magnetic field. Let us consider 
the last open field line of the dipole magnetic field. The photons emitted by backflowing positrons 
produce pairs on adjacent innermost field lines, which may reduce the effective area of the PC 
subject to the heating. Here we shall present very rough estimate of the geometrical filling factor 
for the PC heating qualitatively illustrating this effect. The photons emitted by backflowing 
positron (moving along the last open field line) at height S+ above the effective radius Re 
produces a pair at the effective radius and at the angular distance from the magnetic axis 

r * = (R E + 5+ ) tan 9 - S+ tan A, (62) 

where 6 is the magnetic colatitude of the point at which the photon is emitted, and A is the angle 
between the tangent to the field line at the point of photon emission and the normal to the stellar 
surface (here we may justifiably neglect the surface curvature). In a flat space limit elementary 
geometrical consideration yields the relationship A = fj, + 9, where [i is the angle between the 
tangent to the field line and the radius-vector of the point on the field line with the magnetic 
colatitude 9. For a dipole field, tan ~ ^tan# ~ i#, so that tan A ~ and we get 

r * = flE ( 1 -£) ( '' < 63> 

bmce V = U pc {{RE + S+j/R} 1 / 2 , where 9 pc is the angular size of a standard PC, we can write for 
the effective angular size (at the actual stellar radius) of the heated PC 




(64) 



2Re) 

Thus, in the case where 5+ ~ R, and Re ~ (1 — 2)R, we get 

0; C «(O.7-O.9)0 pc . 

The effective area of the heated PC may therefore be up to 50 % less than the area of a standard 
PC. The above estimate can be additionally justified by the fact that increases toward the PC 
rim, while Re remains practically constant. We thus suggest that for some pulsars this effect may 
be worthy of discussion. 



Self-Limited Acceleration Zone 
3.1. Numerical Results 



Following the procedure outlined in Sections 24-2J3, we have made numerical calculations 
of the PFFs and the total acceleration length, S c , from equation (|l|), for both upward moving 
electrons and downward moving positrons. The parameters of the thermal radiation from the NS 
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surface are somewhat uncertain. Although measured temperatures of pulsed X-rays lie in the 
range 5 x 10 5 — 10 6 K, the distribution of the radiation is not known. For the calculations of 
this Section, we have assumed that the radiation is isotropically emitted from a hot PC of radius 
Rt = WqR and that R = 8 km. Other possibilities, and how they might affect these results will 
be discussed in Section ||. An electron is started at height ho above the NS surface with 7~ = 1. 
For a given value of ho and thus, Re = R + ho, a "first guess" value of h c , and thus also of S~, 
sets the initial acceleration length, equation (§) is integrated in discrete steps upward from the 
starting point, computing En (from either equation jn| or p3fl), \ c and 7 ca at each step. At 
each step, the pair attenuation lengths, S p (e), of both CR and ICS test photons radiated by the 
particle of energy 7(5) are computed from equation (||). The pair attenuation length, and thus the 
value of S~, also computed at every step, is initially infinite, because the energy of the photons 
is small, but decreases with distance as the energy of the radiated photons increases. Although 
the photon attenuation length continues to decrease, the particle acceleration length is increasing 
and S~ has a minimum. This minimum value of h c = ho + (S~) m i n is adopted as the new value 
of the electron PFF for the assumed value of ho- The electron is accelerated again with the new 
value of h c , producing a new PFF at the next value of h c = ho + {S~) m i n . The process is repeated, 
converging to a self-consistent value of h c . 

When ho <C R, we find that the value of S~ for ICS photons is smaller than that for 
CR photons. PFFs at low altitudes are therefore ICS-controlled. Figure § shows examples of 
self-consistent solutions of the ICS-controlled electron PFF near the NS surface (left panel) and the 
corresponding ICS-controlled positron PFF (right panel). In the left panel of Fig. ||[ an electron 
starts very near the NS surface (at height ho = 0.01 R), where the density of soft thermal photons 
and the loss rate due to ICS is high. In this case, photons from resonant ICS produce pairs well 
before CR photons and thus the ICS photons establish the PFF. In fact, the loss rate for CR is 
orders of magnitude smaller than that for ICS at the PFF. In the right panel of Fig. [21 the loss 
rates and integrated energy of the positron, accelerated by the same electric field as the electron, 
is shown as a function of the acceleration length. It this case, however, the larger angles between 
the thermal X-ray photons and the positrons allow them to accelerate through the resonant part 
of the ICS losses at lower energies, and non-resonant ICS of photons above the cyclotron energy 
produce the PFF. Since at pulsar field strengths, scattering above the cyclotron energy occurs in 
the KN regime, the scattered photon energies are much larger than those of photons scattered 
in the resonance. Therefore, the ICS photons producing the positron PFF have higher energies 
than the ICS photons producing the electron PFF (see Fig [?]) and the positron PFF forms in 
a shorter distance. This is the major cause of the difference between the electron and positron 
ICS-controlled PFFs, which will be discussed in more detail in connection with Fig. |8|. Figure ^ 
shows an example of the self-consistent solution for a CR-controlled electron PFF. The electron 
begins accelerating at a significant fraction of a stellar radius above the surface, where the soft 
photon density has dropped, CR losses dominate and CR photons establish the PFF. Because CR 
is much less efficient than ICS, the particles must accelerate to higher energies to produce pairs 
and the total acceleration length S c is longer. 
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We first explore the solutions for the electron PFF, as a function of pulsar parameters, for 
acceleration from the surface. Figures |] and [5] show the height of the electron PFF and the 
maximum acceleration voltage (7 max ) as a function of the scaled colatitude £ for CR-controlled 
PFFs, assuming that ICS is "turned off'. These results can be directly compared with Fig. 5 of 
A83, who computed the PFF due only to curvature photons. Figure |2| shows the strong effect 
of obliquity x on t ne shape of the electron PFF and on the acceleration voltage, due primarily 
to the dependence of En on \- Since the frame-dragging component of En proportional to cosx 
(cf. equation |18|| ) is much stronger near the surface than the component proportional to sinx, 
it dominates at small and intermediate values of x- The frame-dragging electric field becomes 
comparable to that in a flat space for orthogonal rotators (x = 90°), where our solutions match 
very well with those of A83's Fig. 5. The height of the CR-controlled PFF (S~) increases and 
7max decreases with increasing as En decreases. The height of the PFF increases sharply (and 
in fact goes to infinity) both at the magnetic pole (£ = 0) and at the PC rim (£ = 1). These 
features, referred to as slot gaps by Arons, form at the pole and equator due to different effects. 
The gap at the pole occurs because the radius of curvature of the magnetic field lines is infinite 
there, causing the pair attenuation length to go to infinity. The electrons keep accelerating to high 
altitudes and thus reach high energies. At the PC rim, the boundary condition on the potential 
($ = 0), screens En close to the rim, preventing electron acceleration to high enough energy to 
produce any pairs. Thus, 7 max also goes to zero at the rim. Figure || shows the dependence of the 
PFF height and acceleration voltage on surface magnetic field strength B. While the height of the 
PFF decreases with increasing B, due both to an increase in En and a decrease in pair attenuation 
length, 7 m ax is roughly constant even though B varies over two decades. This is because the longer 
total acceleration distance for lower field strengths compensates for the lower rate of acceleration. 

Figure |6| shows examples of ICS-controlled PFFs near the NS surface. We find that both the 
height of the PFF and the acceleration voltage is much lower than for the CR-controlled PFFs. 
This is due to the higher efficiency of ICS in producing photons of pair-producing energy. ZQLH97 
obtained a similar result for ICS-controlled PFFs using the electric field of RS95. Electrons having 
Lorentz factors of only 10 5 — 10 6 , depending on £, are capable of radiating photons of about 
10% to almost 100% of their energy, and these photons will produce pairs. Again, the maximum 
Lorentz factor, 7 maX) is a very weak function of field strength. The ICS-controlled PFF height, 
however, behaves very differently near the magnetic pole. Instead of forming a slot gap, like the 
CR-controlled PFFs, the height of the PFF actually decreases near the pole. ICS photons are 
radiated at much larger angles to the field and thus have shorter pair attenuation lengths than 
curvature photons. 

The above determinations of the electron PFF are not fully self-consistent in that they have 
neglected pairs produced by the positrons that slow down to screen E\\ and accelerate down toward 
the stellar surface. Even if the number of returning positrons is small compared to the number of 
primary electrons, the multiplicity of the downward cascades is high (DH96), due to the increasing 
field strength. From our simulations of downward going cascades using the code of DH96, we find 
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that the multiplicity of a positron accelerated toward the PC from a height of one stellar radius at 
the PC rim is M p ~ 10 4 — 10 5 . If the fraction of positrons that are accelerated downward from the 



electron PFF is f p = A + /A_ ~ 3 x 10 -3 — 2 x 10 _1 (see equations [58] and Q]), then there are 
fpMp ~ 30 — 2 x 10 3 pairs produced by downward cascades for each primary electron. This would 
seem to be sufficient to screen E\\ and form a lower PFF. 

To compute the location of the lower PFF, a test positron is started at the upper (electron) 
PFF with an energy of 7 + = 1. In reality, there is a spread in energies of the pairs created near the 
electron PFF, so that the most energetic ones will decelerate but not turn-around, while the less 
energetic ones will decelerate before they reach the PFF and accelerate downward with varying 
initial energies. We have thus neglected the detailed kinematics of this turn-around process. 
Equation (2) is then integrated to follow the downward acceleration of the positron and its energy 
losses due to both ICS and CR, in the same manner as for the electron. Likewise, the attenuation 
lengths of the photons and the total acceleration length of the positron are computed. Although 
the PFFs of the positrons, like the electrons, are produced by pairs from the ICS process near the 
NS surface in the presence of thermal radiation from a hot PC, we find that the positron PFF 
is always above the start of the electron acceleration. That is, pair front formation by the ICS 
process is not symmetric for upward and downward going particles. This is because the positrons 
scatter the thermal photons at much larger angles than the electrons. While the ICS photons 
producing the electron PFF scatter in the cyclotron resonance to an energy of ~ 7-B', photons 
producing the positron PFF scatter above the cyclotron resonance, in the KN regime (see Fig. 
H), where the scattered energy is ~ 7. This asymmetery will occur for field strengths B'<0.5, 



above which positrons will form PFFs through resonant ICS. As discussed further in Section 3.2 
stable ICS-controlled double PFFs may exist for high-field pulsars because the scattered energy 
for resonant scattering is independent of incident angle. Figure ^ shows the difference between the 
electron, 7~ in , and positron, 7^ in , energies that produce scattered photons of energy e+ in and e min 
that form the minimum PFFs, as a function of ho- Since their ICS is more inefficient, the electrons 
must accelerate to much higher energies in order to produce scattered photons that pair-produce. 
The efficiency of electron scattering, £min/7min' decreases with height, ho, above the surface, 
while the efficiency of positron scattering increases. Therefore, the positrons are able to radiate 
pair-producing photons after traveling a shorter acceleration path. Since the scattered energies of 
the positrons, at a given height, are greater that those of the electrons, these photons have shorter 
pair attenuation lengths. Thus, the total positron acceleration length, S~ = 5'a(7 m i n ) + ^p( e mm) 
is less than the total electron acceleration length, S£. 

Figure [8| shows total positron and electron acceleration lengths, S£ and S~ as a function of 
the electron starting height ho, for both ICS- and CR-controlled PFFs. In these calculations, 
we have allowed the self-consistent PFFs of electrons and positrons to form with only one of the 
radiation mechanisms operating. When only ICS produces pairs, S£ is always significantly less 
than S~ and the difference increases with height. However, when ICS is suppressed and only CR 
is allowed to operate, S£ and S~ are equal. That is because CR for electrons and positrons is 
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the same, the only difference being in radius of curvature between the radiation points, which is 
negligible since S~ <C R. While ICS cannot provide self-consistent solutions to the double PFFs, 
CR can. We suggest that stable, self-consistent double PFFs can only exist when they are formed 
by CR. Since ICS will dominate near the surfaces of NSs radiating thermal soft X-rays, evidence 
of which has been observed in many pulsars, stable PFFs can only form at a height above the 
surface where CR becomes dominant. We take a rough estimate for the height of CR-control of 
the electron and positron PFFs to be where CR energy loss of the electrons dominates over ICS 
energy loss. If electrons of energy 7~ in radiate more curvature photons than ICS photons, then 
the ICS PFF will disappear, and the electrons will continue accelerating until a curvature PFF 
is established. When the positrons accelerating downward from the electron curvature PFF also 
establish curvature photon PFFs, then stable acceleration can occur. In reality, the PFFs will 
switch from ICS to CR control when the number of pairs produced by ICS photons becomes too 
small to screen Eh. But to determine the number of pairs necessary for establishing a lower PFF, 
we would need to know the number of returning positrons. Figure ^ shows an example of the 
electron energy loss rates due to ICS and CR at energy 7~ in as a function of height of the lower 
PFF. At low altitudes, ICS losses dominate by many orders of magnitude, but decrease with ho 
mostly due to the decrease in the thermal photon density (cf. equation [40). The energy 7~ in of 
electrons that are forming the ICS PFF increase with altitude due to the decreasing efficiency of 
ICS. Their CR losses, which are a strong function of energy, therefore increase. The height at 
which electron ICS and CR losses are equal, which roughly set the location of stable acceleration, 
will depend on the PC thermal temperature and radius (cf. equation ]76fl). 



We have computed the height of the stable acceleration zones where CR losses dominate 
control of PFFs for various pulsar parameters. As shown in Figure 1C, the location of the lower 
(positron) PFFs ho are at higher altitudes for higher surface magnetic field strengths, because the 
ICS losses are proportional to B 2 for resonant scattering. The height of the lower PFF decreases 
for increasing pulsar rotation period, P, because at a given colatitude £ the radius of curvature 
increases, but the electron energy 7~ in increases so that CR can dominate at a lower altitude. 
Figure [ll] shows the width of the CR-controlled acceleration zone, S~ (i.e. the distance between 
the lower and upper curvature PFFs) and the maximum voltage drop, 7 max > as a function of pulsar 
period and surface field strength. The maximum particle energy is again remarkably insensitive 
to pulsar parameters, even though the acceleration zone width shows substantial variation. The 
width S~ tends to be smaller for larger fields because both the acceleration length and the pair 
attenuation length are shorter. It has a large increase at long periods, where the pair attenuation 
lengths become large due to the increase in the PC radius of curvature. In fact, there will be a 
maximum period for which electron PFFs can form at a given colatitude. 

Figure [H| shows solutions for the width and maximum Lorentz factor of the self-consistent, 
stable acceleration zones that are controlled by CR double PFFs. The width of these zones, 
located nearly a stellar radius above the surface, are larger than CR-controlled zones at the 
surface (e.g. Fig. || and ||), because the magnetic field has fallen from its surface value, so that 
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E\\ is lower, and the radius of curvature is larger. However, the maximum Lorentz factor of the 
higher altitude acceleration zones is roughly the same as that at the surface. Again, the size of the 
acceleration path has adjusted itself to maintain the same maximum energy, which depends only 
on the geometry of the magnetic field. 



3.2. Analytic Estimates 

It is difficult to derive accurate analytic expressions for the acceleration zone parameters we 
have computed numerically, such as height and width of the stable accleration zone. The form of 
the electric field changes as a function of acceleration length (see Appendix), starting as a linear 
function and then saturating at s ~ -R pc to become a constant. Consequently, there are no simple 
expressions for the processes that apply for altitudes from the surface to several stellar radii. 
Nevertheless, we can characterize the behavior of the solutions in different regions. An estimate for 
ho, the height at which control of the PFF formation switches from ICS to CR, can be obtained by 
finding the height at which the loss rates for CR and ICS of positrons with energy 7 m j n are equal: 

7ic(7min) = 7 CR (7min)- (65) 

We can obtain an estimate of 7m i n by solving equation ([[]) . For s < 3i? pc , we can use the expression 
for Eu given in equation (A3) of the Appendix, which is linear in s = zRe and thus gives a 
quadratic energy increase in particle energy with s. The particle acceleration length is then: 

5 a ( 7 ) = (7/A) 1/2 (66) 

where 



A ^-\2mcV \ c jl-eUJ \Re 



2 

cosx- (67) 



Note that we have rescaled equation (A3) to the effective radius Re, at which the acceleration 
begins. The pair attenuation length can be approximated as the mean free path 

valid for B' < 0.1. Equation (g) may then be written, 

5 c = min{( 7 M 7 )V2 + ^|, ( 69) 

where 



and 



(71) 
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We wish to estimate 7 m i n when ICS produces the PFFs. For nearly all pulsar field strengths 
{B <0.5 B cr ) the photons which produce the positron PFFs have been blue-shifted above the 
resonance and scatter in the KN regime (see Fig. |2|). In this case, we can approximate the 
scattered photon energy as e ~ 7. Inserting this value of e into equation (|69|), we then obtain 7min 
giving the minimum S c by setting the derivative of S c with respect to 7 to zero: 

72 =(2^-) =3.4xlO»B^(^)(o„s,)V3. (72) 

where B\2 is the surface field strength. Thus, the particle energy required to radiate pair-producing 
photons increases with altitude, as was found in our numerical results (cf. Fig. [?]). 

We can obtain an estimate for the width of the ICS controlled acceleration zones by 
substituting 7^^ from equation ([72]) into equation (^) for the minimum value of S c , 

S? N ~ 5.3 x 10 3 B u 2/3 P 1 ' 2 R\ 12 (^f (cos x)' 1/3 cm. (73) 

We can estimate the height ho above the surface at which CR energy loss exceeds ICS losses for 
the downward moving positrons of energy 7 m i n . From our numerical results, we found that at 
the altitudes where CR begins to dominate the particle energy loss, the ICS is in the K-N regime 
above the resonance. We will therefore equate the curvature loss rate to the K-N loss rate to 
obtain an estimate for ho- Substituting 7 m ; n from equation ([T2]) in equation (^) for the CR energy 
loss, we have 

TaaCff) = 1.6 x ltfB^R-^P^ (^) 3 (cos^s" 1 (74) 

which sharply increases with altitude, Re = R + ho. We approximate the K-N loss rate, modifying 
Blumenthal & Gould's (1970) formula for scattering of isotropic thermal photons of temperature 
Tq, from a hot polar cap of radius Rt, with a dilution factor that accounts for the fall-off in photon 
density with height h: 

.JCN\ _ ^0 (mckT) 2 ^ ( 4-f min kT\ R\ _ _ n R ^ 1 n8 ^ 2 R T 



7kn (7min ) - 7^5 7 In ' ' ,0 9 - 7 - 6 X 10 S ~ T 6 TT (75) 

/KNV/mm7 gmc 2 ^3 \ mc 2 J h 2 + Rl 6 h 2 V ' 



for Rt <C h and ignoring the slowly varying In factor. Setting 7 CR (7 m ^^ r ) = T^C^mta) an d 
assuming Re ~ ho, 

|) ^ - 8.6 BT P 1/5 RT tT (f ) (cos x)- 4/15 - (76) 

Although the above expression involves a number of approximations, it does seem to roughly 
reproduce our numerical results. For example, taking the parameters Tq = 0.5, Rq = 0.8, P = 0.1, 



B\2 = 4.4 and Rt = 0.04 R gives ho/R — 1.6, fairly close to the value shown in Fig. 11 



At the altitude ho, where CR losses are dominant and curvature photons produce the PFFs, 
the width of the acceleration zone can be calculated in a similar manner to that of the ICS 
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controlled acceleration zones. Using the CR critical frequency, e cr = (3/2)(h/mc)'y 3 / p c in equation 
(p9|), we can solve for the particle energy 7^5 which produces the minimum acceleration zone 
width S c , 

tS£ = (6C P < 2 ) 2 / 7 = 4.7 x itfB^pW (§) 4/7 ^ /? (cosx) 1/7 , (77) 



where 

r..-( 

Then substituting the expression for back into equation (6£) for S c , we have 



C P =(^). (78) 



S C C R ~ 4.8 x 10 4 B x - 2 4 / 7 P 4 / 7 < 7 (^f ) 16/7 (cos X )- 3/7 cm. 



(79) 



The expression for 7^j? above is also a very good estimate of 7^, the maximum Lorentz factor 
of the CR-controlled accleration zone, because S p <C S a for CR photons (i.e. the PFF is very close 
to where the first pairs are produced). 

For very high field strengths (B' ~ 0.5) both electron and positron ICS operate in the resonant 
scattering regime. The incident soft photons scatter in the cyclotron resonance and their scattered 
energies are approximately e s ~ jB' . In the Thomson limit of resonant scattering, the energy of 
the scattered photons is independent of incident photon angle. This means that = S~, PFF + 
will coincide with ho, and stable ICS-controlled acceleration zones are thus possible for high-field 
pulsars. However, in the relatively high fields where resonant ICS controls the positron PFFs, 
relativistic effects on the resonant scattering cross section as well as photon splitting will become 
important (as will be discussed in Section ||) and our calculation of PFFs is incomplete. 



4. Discussion 

In this paper, we have investigated the effect of cascades from downward-accelerated positrons 
on the electrodynamics of the PC particle acceleration. We find that when ICS produces pairs in 
the acceleration zone, the positron cascades may screen the accelerating electric field and disrupt 
particle acceleration near the NS surface. Thus, if lower PFFs can develop, the picture of steady 
particle acceleration from the PC surface must undergo major revision. We suggest that a stable 
acceleration zone may exist at an altitude of about one stellar radius above the PC, in which pairs 
from CR limit the electrostatic acceleration of primaries by screening the electric field near the 
upper PFF. 

The calculations presented here are only a first attempt to describe the physics of what is a 
very complicated process. We have made many assumptions and approximations to obtain our 
results. While we believe that the gross qualitative results of our study are correct, there are a 
number of aspects which should be treated more accurately to achieve more solid quantitative 
results. We have assumed that the screening of the Eh occurs over the short distance determined 
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by the upper boundary condition of Poisson's equation, not by the dynamics of the pair screening. 
The details of the pair screening also determine the fraction of pairs that return to the PC and we 
have briefly outlined in Section ( |2.4[ ) how such a calculation can be done within the electrodynamic 
framework set up in this paper. A determination of the returning positron fraction would answer 
some interesting questions, such as whether PC heating is important relative to cooling in setting 
the PC temperature and if so, whether there is a feedback loop between the formation of the 
upper PFF and the PC heating. The characteristics of the downward cascades need further study 
and modeling. While we have done preliminary modeling of these cascades near the PC rim, there 
will be significant variation in pair yields with magnetic colatitude. Near the magnetic pole, the 
pair multiplicity will drop, allowing the lower PFF height to decrease or even disappear. 

We have found that the nature of the PC acceleration depends strongly on the characteristics 
of the radiation from the hot PC. Thus, to solidify the quantitative aspects of our results, it is 
important to treat the thermal PC radiation and the ICS process as accurately as possible. For 
example, we have assumed in this paper, that the thermal radiation is uniformly emitted over a 
PC of size Rt = 3#o with an isotropic flux distribution. However, studies of thermal radiation 
propagating through a strongly magnetized NS atmosphere (Pavlov et al. 1994) will not be 
isotropic, due to the anisotropy of the magnetized scattering cross section. The expected radiation 
pattern consists of a pencil component, beamed along the magnetic field and a fan component 
perpendicular to the field. Such a beam pattern has been found to be consistent with observed 
thermal X-rays pulse fractions and pulse profiles for several pulsars (Shibanov et al. 1995, HM98). 
In addition, our treatment of ICS uses a combination of the (non-relativistic) magnetized cross 
section in the Thomson limit for scattering near the fundamental cyclotron resonance, and the 
(non-magnetic) Klein-Nishina cross section to describe relativistic effects above the resonance. 
While this hybrid treatment is somewhat inaccurate for B ~ 0.1£> cr and therefore not completely 
satisfactory, the fully relativistic QED magnetic scattering cross section (e.g. DH86) is too 
complicated for use in this type of calculation. In particular, the quantization of the electron 
momentum perpendicular to the field limits the number of Landau states contributing to the cross 
section for each incident photon energy. The magnetic QED cross section for scattering just above 
the fundamental cyclotron resonance, and therefore the loss rate, could therefore be substantially 
lower than the Klein-Nishina cross section we have used here. Unfortunately, there exists no 
simplified, approximate expression for the QED scattering cross section, which smoothly bridges, 
and allows a unified treatment of, the relativistic resonant and non-resonant regimes of ICS. 

A physical process which we have neglected in this study, photon splitting, is not expected to 
be significant for the magnetic fields we have considered, but will be very important for pulsars 
having surface Bq ~ 0.5 -B cr . Photon splitting, a third-order QED process in which one photon 
splits into two, operates only in very high magnetic fields and competes as an attenuation process 
with one-photon pair production (Harding et al. 1997) because it can occur below pair threshold. 
The implications of photon splitting for PC PFFs is profound. Pulsars having Bq ~ 0.5 B cr will 
produce fewer pairs, especially near the surface, so that the cascades from downward- moving 
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positrons may not produce a lower PFF. This is about the same field strength where we found 
that double PFFs controlled by resonant ICS at the NS surface become possible. Thus, the 
stable, CR-controlled double PFF structure we have studied in this paper, which we found to 
move to higher altitudes with increasing surface field strength, will eventually collapse back to the 
surface at very high B$. A single PFF will then form, controlled by ICS, and the acceleration 
zone will have characteristics similar to that shown in Figure H At extremely high surface fields, 
Bq ~ 1.0 -B CT photon splitting will suppress pair creation completely at this single PFF, and 
the pulsar may be radio quiet (Baring & Harding 1998). However, the particle acceleration in 
these pulsars will operate very efficiently, free of any screening of Eu, so they are expected to be 
observable at high energies. Bound-state pair production, where photons convert to positronium 
just below pair threshold rather than converting to free pairs (Usov & Melrose 1995), may also 
come into play at higher field strengths. 

Some of the main results that we have presented in this paper have some important 
implications for pulsar high-energy emission. One of these results is the insensitivity of the 
acceleration voltage (maximum particle energy) 7 max to any pulsar parameters such as period, 
surface magnetic field strength, obliquity (except for nearly orthogonal rotators), and even height. 
This result showed up many times in the course of our calculations (cf. Figs. f|, |5j, || [ll] and 
12, and equation j77|), and seems to be a robust characteristic of this type of PC acceleration 



model. The maximum particle acceleration energy varies only within each pulsar, as a function 
of magnetic colatitude. This energy, between 5 x 10 12 eV and 5 x 10 13 eV for CR-controlled 
acceleration zones, is about two to three times higher than that without frame-dragging (e.g. A83), 
and is consistent with the primary particle energy required in CR-initiated PC cascade models of 
7-ray pulsars (e.g. DH96). The acceleration energy is very high near the magnetic poles, where 
the radius of curvature goes to infinity, allowing for the possibility of a narrowly beamed, hard 
component in 7-ray pulses. However, the radiation power emitted by these high-energy particles 
will be small, because the curvature radiation loss rate is proportional to p~ 2 . The insensitivity 
of 

7max to pulsar parameters implies that the primary particles in all pulsars are accelerated to 
the same energy, and that the luminosity of the high-energy emission should depend only of the 
flux of primary particles. This is consistent with trends in the observed 7-ray pulsar luminosities 
(Thompson et al. 1997). The 7-rays will originate within a stellar radius of the upper stable PFF, 
at height h c = ho + S~ . From our results of Figs. [U^ and |ll|, and equations (^) and (ff|), the 
height of the 7-ray production increases with period, roughly as P 1 / 2 . The standard PC half-angle 
at height P 7 = R + h c will be 6 C ~ (QR^/c) 1 / 2 oc P~ 1//4 . Thus, the 7-ray emission solid angle fly, 
which is expected to be, $7 7 ~ 27r[l — cos(3# c /2)], will be very weakly dependent on period and 
field strength. 

Our conclusion that stable acceleration may occur in most pulsars at some altitude above the 
surface will have consequences not only for high-energy emission, but for radio emission as well. If 
electron-positron pairs are necessary for coherent radio emission, then the dependence of the PFF 
altitude on pulsar parameters should be taken into account when determining the radio pulsar 
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"death line" , the line on the period-period derivative diagram beyond which pulsars are incapable 
of producing pairs. Our calulcations in this paper suggest that pulsars with long periods do not 
produce PFFs, and that there will be a "death line" at periods below where there are observed 
radio pulsars. But this is a long-standing problem of PC acceleration, most recently discussed by 
Arons (1998), and one that should be addressed in future studies. 

Our principal findings can be summarized as follows. 

1. Lower PFFs may form by positrons returning to the NS surface from the upper PFF. 

2. Pair creation by the ICS process dominates near the stellar surface, but is not symmetric for 
upward and downward going particles, so that stable, double PFF formation is very unlikely. 

3. Stable, self-consistent double PFFs can only exist when they are formed by CR. They can 
only form at a height above the surface where CR becomes dominant. 

4. The maximum particle energy is insensitive to any pulsar parameters such as period, surface 
magnetic field strength, obliquity, even height, and is sufficient to power 7-ray pulsars. 

The main conclusion of this paper is that the cascades from positrons returning to the PC 
may have a significant effect on the primary particle acceleration in pulsars and should not be 
neglected. It is possible that lower PFFs do not form for all pulsars, and may not form over the 
entire PC. A detailed study of the screening of the accelerating electric field by the returning 
positron cascades is beyond the scope of this paper, but will ultimately be necessary to understand 
PC acceleration. These studies are needed to address the questions of the returning positron 
fraction and the multiplicity of downward cascades. If we can show that the returning positrons do 
not screen E\\, then the present assumption of acceleration right from the NS surface is valid. But 
if these studies show that screening at a lower PFF is effective, then the possibility of acceleration 
above the NS surface must be incorporated in PC models. 

We thank the referee Bing Zhang for his very careful review and insightful comments. We are 
also grateful to Joe Daugherty, for help in the downward cascade simulations, to Steve Sturner for 
discussions on inverse-Compton scattering and Matthew Baring for comments on the manuscript. 
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APPENDIX A. APPROXIMATE EXPRESSIONS FOR THE ACCELERATING 

ELECTRIC FIELD 

Here we present the explicit expressions for the electric field component parallel to the 
magnetic field for the various ranges of altitude. These expressions have not been rescaled to allow 
for acceleration starting at an effective radius above the stellar surface, so that z is the altitude 
above the surface. 



S c < r pc , and 0.1 r pc (r pc /R) < s < S c : 
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In the above expressions r pc ~ Qq R is the PC radius, z c = S c /R is the dimensionless altitude of 
the upper boundary (PFF), and all other quantities are scaled with the true stellar radius. Note 
also that fl"(l)<S(l) ~ 1. 



-33- 



REFERENCES 

Arons, J. 1981, in IAU Symp. 95, Pulsars, ed. W. Sieber k Wielebinski (Dordrecht: Reidel), 69 
Arons, J. 1983, ApJ, 266, 215: A83 

Arons, J. 1998, in Neutron Stars and Pulsars, eds. N. Shibazaki, N. Kawai, S. Shibata k T. Kifune 
(University Academy Press: Tokyo), p. 339. 

Arons, J., k Scharlemann, E. T. 1979, ApJ, 231, 854: AS79 

Baring, M. G. k Harding, A. K. 1998, in preparation. 

Becker, W., k Triimper, J. 1997, A&A, 326, 682 

Bussard, R. W., Alexander, S. B., k Meszaros, P. 1986, Phys. Rev. D, 34, 440 
Canuto, V., Lodenquai, J., k Ruderman, M. 1971, Phys. Rev. D, 3, 2303 
Daugherty, J. K., k Harding, A. K. 1983, ApJ, 273, 761: DH83 
Daugherty, J. K., k Harding, A. K. 1986, ApJ, 309, 362: DH86 
Daugherty, J. K., k Harding, A. K. 1989, ApJ, 336, 861. 
Daugherty, J. K., & Harding, A. K. 1996, ApJ, 458, 278: DH96 
Dermer, C. 1990, ApJ, 360, 197: D90 
Deutsch, A. 1955, Ann. d'Ap., 18, 1 

Fawley, W. M., Arons, J., k Scharlemann, E. T. 1977, ApJ, 164, 529 

Goldreich, P., k Julian, W. H., 1969, ApJ, 157, 869 

Gonthier, P. L. k Harding, A. K., 1994, ApJ, 425, 767 

Harding, A. K., Baring, M. G. k Gonthier, P. G. 1997, ApJ, 476, 246. 

Harding, A. K., k Muslimov, A. G. 1998, in Neutron Stars and Pulsars, eds. N. Shibazaki, N. Kawai, 
S. Shibata k T. Kifune (University Academy Press: Tokyo), p. 311. ( |astro-ph /9802038 ) 

Harding, A. K., k Muslimov, A. G. 1998, ApJ, in press. (|astro-ph/9802044|) 

Jones, P. B. 1978, MNRAS, 184, 807 

Jones, P. B. 1985, Phys. Rev. Lett., 55, 1338 

Jones, P. B. 1986, MNRAS, 218, 477 



-34- 



Kardashev, N. S., Mitrofanov, I. G., & Novikov, I. D. 1984, Sov. Astron., 28, 651 
Lorentz, C. P., Ravenhall, D. G., & Pethick, C. J. 1993, Phys. Rev. Lett., 70, 379 
Mestel, L. 1998, Stellar Magnetism (Oxford: Oxford University Press) 
Muslimov, A. G., & Harding, A. K. 1997, ApJ, 485, 735: MH97 
Muslimov, A. G., & Tsygan, A. I. 1990, AZh, 67, 263: MT90 
Muslimov, A. G., & Tsygan, A. I. 1992, MNRAS, 255, 61: MT92 
Neuhauser, D., Langanke, K., & Koonin, S. E. 1986, Phys. Rev., A33, 2084 
Neuhauser, D., Koonin, S. E., & Langanke, K. 1987, Phys. Rev., A36, 4163. 
Pavlov, G. G., Shibanov, Yu A., Ventura, J. & Zavlin, V. E. 1974, A & A, 289, 837. 
Ravenhall, D. G., & Pethick, C. J. 1994, ApJ, 424, 846 
Ruderman, M. A., & Sutherland, P. G. 1975, ApJ, 196, 51: RS75 

Shibanov, Yu. A., Zavlin, V. E., PAvlov, G. G. & Ventura, J. 1995, in 17th Texas Symposium on 
relativistic Astrophysics, eds. H. Bohringer, G. Morfill & J. Trumper, Annals of NY Acad. 
Sci., 759, p. 291. 

Sturner, S. J. 1995, ApJ, 446, 292: S95 

Sturner, S. J. & Dermer, C. D. 1994, ApJ, 420, L79. 

Sturrock, P. A. 1971, ApJ, 164, 529 

Thompson, D. J., Harding, A.K., Hersen, W. & Ulmer, M.P. 1997, in Proc. of the 4th Compton 
Symposium, ed. CD. Dermer, M.S. Strickman & J.D. Kurfess (AIP 410: New York), 39. 

Tsai, W.-Y. & Erber, T. 1974, Phys. Rev. D, 10, 492. 

Usov, V.V. & Melrose, D. B. 1995, Aust. J. Phys., 48, 571. 

Wang, F. Y.-H., Ruderman, M., Halpern, J. P., & Zhu, T. 1998, ApJ, in press 

Xia, X. Y. et al. 1985, A&A, 152, 93. 

Zavlin, S.E., Shibanov, Yu. A. k Pavlov, G.G. 1995, Astron. Lett., 21, 149. 
Zhang, B., & Qiao, G. J. 1996, A&A, 310, 135 

Zhang, B., Qiao, G. J., Liu, W. P., & Han, J. L. 1997, ApJ, 478, 313: ZQLH97 



This preprint was prepared with the AAS IATgX macros v4.0. 




Fig. 1. — Schematic illustration of the double pair formation fronts (PFF) produced by electrons 
and positrons at height ho above a pulsar PC. Parallel electric field accelerates particles in the 
shaded zone, up to the electron PFF - at height h c , and is screened (E • B = 0) everywhere else. 
Upward accelerating electrons of energy 7~ in radiate photons of energy g~ in that can produce the 
first electron-positron pairs to form the PFF - . Downward accelerating positrons of energy 7^ in 
radiate photons of energy e~ in that produce pairs which may screen the electric field at PFF + . 
This configuration shows double ICS-controlled PFFs, which are not stable because the positron 
PFF+ lies about h . 
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Fig. 2. — Parallel electric field En , particle energy 7 and energy loss rates 7 due to non-resonant and 
resonant ICS and CR as a function of acceleration length s for ICS-controlled PFFs of electrons (left 
panel) and positrons (right panel). Here, Kq is the height of the lower PFF (start of acceleration) 
above the surface, Tq = 0.5, P = 0.1 s, B = 0.1 B cr , \ = 0-2 (radians), and £ = 0.7. 
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Fig. 3.- 



Same as Fig. [2], but for the electron CR-controlled PFF at height ho = 0.5R. 
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B = 0.1B cr , P = 0.1, h = 




Fig. 4. — Profiles of the acceleration voltage (maximum electron Lorentz factor) 7 max and width, 
S 1 ", of the acceleration zone formed by a CR-controlled PFF at the NS surface, neglecting losses 
and pairs from ICS, as a function of magnetic colatitude scaled to the PC half angle £ = 6/8q, for 
different obliquities %■ 
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Fig. 5. — Same as Fig. || but for different values of the surface magnetic field strength. In Figs. 
5-10, x is i n radians. 




Fig. 6.— 



Same as Fig. |5| but for ICS-controlled PFFs at the NS surface. 
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Fig. 7. — Minimum Lorentz factors of electrons, 7~ in , and positrons, 7^ that radiate pair- 
producing ICS photons of energy e~ in and e^ in , as a function of acceleration starting height ho. 
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Fig. 8. — Total positron (5+) and electron (S c ) acceleration lengths for CR-controlled and ICS- 
controlled PFFs as a function of acceleration starting height /iq- 
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B = 0.1 B cr , P = 0.1 s £ = 0.7 % = 0.2 




-|0" 2 2 3 4 5 -|0" 1 2 3 4 5 -j Q° 



h /R 

Fig. 9. — Electron energy loss rates for CR and ICS at Lorentz factors 7 min , as a function of 
acceleration starting height ho- 
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Fig. 10. — Height at which electron energy losses, at Lorentz factors 7^i n , from CR and ICS are 
equal, above which the electron PFF is controlled by CR, as a function of pulsar period and surface 
value of magnetic field strength. 
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Fig. 11. — Acceleration voltage (maximum electron Lorentz factor) 7 max and width, S~ , of the 
acceleration zone formed by a CR-controlled electron PFF at the height ho where control of the 
electron PFF switches from ICS to CR (see Fig. 10), as a function of pulsar period and surface 
value of magnetic field strength. 
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Fig. 12. — Profiles of the acceleration voltage (maximum electron Lorentz factor) 7 max an d width, 
S~ , of the acceleration zone formed by a CR-controlled PFF at the height ho where control of the 
electron PFF switches from ICS to CR (see Fig. as a function of magnetic colatitude scaled 
to the polar cap half angle £ = 8/6q, for different values of surface magnetic field strength. 



